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Abstract 

Recent x-ray Compton scattering experiments in ice have provided use- 
ful information about the quantum nature of the interaction between H2O 
monomers. The hydrogen bond is characterized by a certain amount of charge 
transfer which could be determined in a Compton experiment. We use ab- 
initio simulations to investigate the hydrogen bond in H2O structures by cal- 
culating the Compton profile and related quantities in three different systems, 
namely the water dimer, a cluster containing 12 water molecules and the ice 
crystal. We show how to extract estimates of the charge transfer from the 
Compton profiles. 
PACS numbers: 78.70.Ck, 71.15.-m, 31.15.Ar 

I. INTRODUCTION 

The hydrogen bond is one of the least well-understood components in the energy de- 
composition that is used to predict the folding of biological complexes such as proteins. 
Its importance stems from its directionality and modest bonding energies midway between 
strong covalent and weak Van der Waals bonds. For this reason the inter-molecular inter- 
action is difficult to characterize. Recently, X-ray Compton scattering measurements in ice 
revealed subtle oscillations in the profile anisotropy, which can be interpreted as reflecting 
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various quantum mechanical aspects of the hydrogen bond |iy2]||[|,|5) . In simple descriptions 
of the hydrogen bond the highest occupied molecular orbital (HOMO) of a proton 

acceptor momomer can be viewed as a lone pair p state of the oxygen, while the lowest 
unoccupied molecular orbital (LUMO) of the donor neighbor molecule is the antibonding 
OH orbital. Therefore, the problem can be reduced to an effective 2 by 2 eigenvalue calcu- 
lation where the HOMO-LUMO mixing can in principle be measured by an angle 9. In this 
effective model the amount of charge transferred can be defined as 

AQ = J | cos6» < r | > +sin6» < r | 1 > - < r | > | 2 d 3 r (1) 

where < r | > is the wavefunction of the lone pair and < r | 1 > is the wavefunction of the 
LUMO at the donor neighbor molecule. When the mixing is small, one obtains AQ « sin 2 9. 
Besides, in this limit, sin9 turns to be approximately proportional to the energy gain due 
to the hydrogen bond. 

Compton scattering gives us information on the (ground state) electronic momentum dis- 
tribution of the system under study. Of particular interest is the width Ap of the momentum 
density over which the valence momentum density falls to zero around a cutoff momentum 
(Fermi momentum in a metal). Since the LUMO contains higher momentum components 
than the HOMO, the variation of the width Ap is clearly related to the mixing 9. Moreover, 
Friedel and Peter [[| provided a useful discussion on the impact in momentum density of the 
mixing of occupied with virtual orbitals. Their arguments can be extended to the present 
case in order to estimate the charge transfer mixing angle. We limit our discussion here to 
the more general concept of charge transfer without invoking the term covalency ||. As 
a matter of fact ionic compounds are also characterized by a charge transfer but they are 
certainly not described as covalent systems. 

This paper is organized as follows. In Sec. II we present the computational method and 
recall the definitions of the quantities studied. Sect. Ill is devoted to estimate the charge 
transfer and to describe its properties in position space. Finally, Sec IV summarizes our 
main conclusions. 
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II. METHOD 



We start with the one-particle density matrix, p(r, r'), which is defined in terms of the 
normalized N-particle wavefunction, as 

p(r,r') = N Jd^*(r,0*(r',0 ■ (2) 

where the integral extends over the coordinates of all other particles. If the many body 
wavefunction is represented by a single determinant, which is true in the case of Hartree 
Fock then the density matrix is idempotent (p = p 2 ) and reduces to a summation over the 
occupied spin-dependent orbitals ipi, i.e. 

N 

p(r,r')=E^(r)^(r'). (3) 

i=i 

The Electron Momentum Density (EMD) is defined as the momentum transform of the 
density matrix, 

/°(P) = / / d 3 rtf¥p(r, r') exp(-ip ■ (r - r')) . (4) 

Note that, in general, the EMD involves off-diagonal elements of the real space density 
matrix. In connection with the Compton spectra, it should be noted that the Compton 
profile (CP), J(p z ), represents the double integral of the ground-state EMD pip): 

J(Pz) = J J piv)dp x dp y} (5) 

where p z lies along the scattering vector of the x-rays. 

The program used in the present work was CRYSTAL98 ||10|| , which is especially appro- 
priate for first-principles calculations of the density matrix p(r, r') in extended (periodic) 
and in molecular systems. For the ice crystal we have used the experimental atomic dis- 
tances. We have also considered the dimer and cluster containing 12 water molecules which 
are shown in Fig. [I]. The 0-0 distances and the orientations of the atoms in the cluster and 
in the dimer are equal to the corresponding ones in the crystal. The CRYSTAL98 program 
employs linear combination of atomic orbitals for solving the Hartree-Fock equation. One 
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can show that the wavef unction involved in the hydrogen bond are described reasonably 
well already at the level of the Hartree-Fock approximation |jn|JI2| . The occupied orbitals 
used to determine the momentum density and the Compton profile are calculated using the 
Restricted Hartree-Fock (RHF) scheme with the 6-31G** basis set ||12|| . 



III. RESULTS 

The density of state (DOS) of the ice crystal is shown in Fig.^|. The band gap is about 
15 eV. In the occupied bands, one can notice two high peaks near -20 eV and near -13 eV. 
For comparison, in the dimer the calculated RHF energies of the 8 valence occupied states 
are -37 eV, -36 eV, -20 eV, -19 eV, -16 eV, -15 eV, -14 eV and -13 eV. The electronic density 
of these dimer states has been visualized by Chaplin [Tj|. One can see that the densities of 



the orbitals corresponding to -19 eV and to -15 eV show overlap across the hydrogen bond. 
Interestingly, the two peaks of the crystal DOS correspond to dimer orbitals orthogonal to the 
hydrogen bond: —20 eV corresponds to a state strongly localized on the acceptor molecule 
while —13 eV corresponds to a state localized on the donor water molecule. Clearly the 
dimer has some very important things to tell us about the nature of the hydrogen bond. 
The projected DOS shown in Fig.[3] reveals that the peaks near —20 eV is a combination 
of 2pO and (1sH,2sH) while the peak near —13 eV is almost a pure 2pO. Moreover, one 
can notice that the highest occupied bands are mostly a combination of 2pO and (1sH,2sH) 
states with a very small 2sO contribution. This picture is therefore consistent with the 
charge transfer theories |7||| explaining the hydrogen bond by the mixing of the oxygen lone 
pair with an antibonding OH orbital centered on a neighbor molecule. 

In order to extract molecular energies and other properties, it can be useful to build a 
density matrix as a superposition of density matrices of isolated molecules arranged in the 
same geometry as in the crystal or the cluster. The keyword MOLSPLIT (in the CRYSTAL 
98 input) performs an expansion of the lattice, in such a way that the molecules of the 
system are at an infinite distance from each other to avoid inter-molecular interactions. Our 
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calculations indicate that the gain in energy per hydrogen bond in the crystal is —0.5 eV. 
This energy gain is strengthened by a factor of two as one goes up in the cluster size. 

Our aim in this work is to probe theoretically the effects of the hydrogen bond in water 
by evaluating the changes it provokes in the electronic momentum distribution which can 
be probed using an experimentally accessible quantity, the Compton profile. For this we 
use the spherical average of the momentum density n(p) which can be obtained from the 
isotropic, or spherically averaged Compton profile J(p) by the formula 

( \ ldJ (p) / R N 
p dp 

The spherical averages J(p) and n(p) are thus used for the evaluation of a quantity Ap, 
defined as the full-width at half maximum of the broad maximum exhibited by the function 
—dn(p)/dp (see Fig. f|). This maximum corresponds to the typical valence cutoff momentum 
(the Fermi momentum in the case of a metal). In the limit of a small mixing angle, formulae 
similar to those given by Friedel and Peter |J indicate that sin 9 is proportional to the 
variation of the momentum density smearing width. Moreover, by assuming a linear scaling 
between the momentum width and the energy separation of the effective 2 by 2 eigenvalue 
problem, we obtain 

1 Ap - Ap 

smt) & —= , (7) 

y/2 Ap ' K ' 

where Apo corresponds to the single water molecule and Ap to hydrogen bonded water 
structure being probed. 

The radial derivative of —dn(p)/dp is shown in Fig. |5| for the water molecule, the dimer, 
the cluster containing 12 water molecules and the ice crystal. The increase of Ap as a 
function of the cluster size is striking. The sin8 given by Eq. (|7]) it is strengthened by a 
factor of two as one goes from the dimer to the crystal. Thus, these estimates yield values 
consistent with the trends observed in the total energy calculations. The amount of charge 
transferred is about 0.5% of an electron in the dimer , 1% in the cluster and 2% in the 
crystal. Interestingly these values are in reasonable agreement with other estimates of the 



charge transfer given in previous quantum mechanical calculations |7]|||T4[]. Finally, in Fig. 
|5] we also notice a shift of the maximum of the function —dn(p)/dp as a function of the 
cluster size, which can explain an observed increase of the kinetic energy. 

In order to show how Compton scattering data can assist in the interpretation of the 
bonding properties in position space, one can define a function B(r) which is the Fourier 
transform of the momentum density |l5|Jl6| , 



B(r) = J p(p) exp(-ip • r) d 3 p . (8) 

From the convolution theorem, B(r) is just the autocorrelation of the one-electron wave 
functions. Moreover B(r) in a given direction can be calculated as the Fourier transform 
of the Compton profile in this direction. The orbital contribution of the Is core electron 
of oxygen vanishes beyond 1.5 A so that the B(r) is dominated by valence electrons. By 
analyzing B(r) in the direction z parallel to the hydrogen bond, we observe a negative local 
minima near 2 A, a positive local maximum near 3 A and finally a negative local minima at 
higher distances as shown in Fig|| The corresponding lines along x and y, also plotted in 
in Fig|], differ substantially from the z direction. Both x and y curves are quite similar and 
present only the local negative minima near 2 A, which is broader toward the large distances, 
and then their amplitudes go to zero at larger distances. The total autocorrelation function 
B(r) can be also interpreted in relation to the hydrogen bond with a comparison to isolated 
aligned molecules (obtained by MOLSPLIT). The isolated water molecules curves are shown 
in the bottom part of Fig.|| along x, y and z. One can observe that none of these curves 
presents the oscillations that characterizes the z direction in the crystal. The autocorrelation 
function B(r) along the direction of the strongest hydrogen bond is shown in Fig. [7|for the 
crystal and the dimer. Both curves deviate strongly from the monomer curve. They also 
present some clear differences among them: the crystal possesses a deeper minimum and a 
somewhat lower maximum at 3 A. However, beyond 4 A the two curves cannot anymore be 
distinguished indicating that the charge transfer occurs mostly within this distance range. 
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IV. CONCLUSION 



In conclusion, this study confirms that the Compton profile is a useful tool for the 
investigation of the quantum nature of hydrogen bond in ice |l| and probably in other 
compounds such as urea \T7 , RNA and DNA [HJ. Our results show a unique sensitivity 
of Compton scattering to the momentum width Ap. In this work an attempt has been 
done to give an estimate of the HOMO-LUMO mixing from Compton profiles. The present 
method can be also used for amorphous materials where Compton profiles anisotropics are 
not available. Moreover, we have used the directional Compton profiles to study the bonding 
properties in the position space by using the autocorrelation functions. We have determined 
that the charge transfer occurs mostly within a range of about 4 A . The present work also 
shows that x-ray inelastic scattering experiments can shed light on the important role in short 
hydrogen bonds played by quantum mechanical charge transfer between nearest molecules. 
These phenomena cannot be described by the standard effective classical potentials used for 



instance in some water clusters simulations 19 
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FIGURES 




FIG. 1. Water cluster containing 12 water molecules. 
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FIG. 2. DOS of the ice crystal according to the RHF calculation. 
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FIG. 3. Projected DOS per orbital of the ice crystal according to the RHF calculation. 
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FIG. 4. Spherically averaged Compton profile J(p) of water (full line) together with the function 



—dn(p)/dp (dashed line) used to define the quantity Ap. 
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FIG. 5. Negative radial derivative of n(p) for the water molecule, the dimer, the cluster con- 
taining 12 water molecules and the ice crystal. 
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FIG. 6. Top: theoretical B(r) for the Ice crystal along x (dashed-dot line), y (dashed line) 
and z (full line) directions obtained from the Compton profiles calculated by CRYSTAL98 fi~(| . 
Bottom: corresponding lines for the isolated molecules. The lines are normalized so that B(0) = 1. 
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FIG. 7. Top: theoretical B(r) for the Ice crystal (full line) for the dimer (dashed line) and the 
water molecule (dashed-dot line) in the direction of the strongest hydrogen bond. The lines are 
normalized so that B(0) = 1. 
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